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Abstract 

We discuss the semiclassical approaches for describing systems with spin-orbit interactions by Little- 
john and Flynn (1991, 1992), Frisk and Guhr (1993), and by Bolte and Keppeler (1998, 1999). We 
use these methods to derive trace formulae for several two- and three-dimensional model systems, 
and exhibit their successes and limitations. We discuss, in particular, also the mode conversion 
problem that arises in the strong-coupling limit. 
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1 Introduction 

The periodic orbit theory (POT) initiated by M. Gutzwiller over three decades ago [|l| has proven 
to be a successful tool for the semiclassical description of chaotic systems [||, ^ |], ||. Several 
extensions of Gutzwiller's semiclassical trace formula to systems with regular and mixed dynamics 
1^, 0, ||, 1C, 11, 12, |l^, [TJ] have made it possible to describe quantum shell effects in many physical 
systems in terms of the shortest classical periodic orbits (see 16] for recent surveys). 

However, in all the approaches mentioned so far, the spin degrees of freedom have not been 
incorporated in the semiclassical theories. This becomes, in particular, necessary when one wants 
to apply the POT to systems with spin-orbit interactions, such as nuclei, atoms, or semiconductor 
nanostructures. Littlejohn and Flynn [jlTj have developed a semiclassical theory of systems with 
multi-component wavefunctions and applied it |l8| to the WKB quantization of integrable spherical 
systems with the standard spin-orbit interaction; Frisk and Guhr [[ll| have extended their method 
to deformed systems with spin-orbit interaction. None of these authors have, however, developed 
an explicit trace formula. Bolte and Keppeler |20[ | have recently derived a relativistic trace formula 
from the Dirac equation. They studied several non-relativistic limits and rederived the Littlejohn- 
Flynn (LF) approach in the limit of a strong spin-orbit coupling, thereby justifying some ad hoc 
assumptions made in [jl9| ]. A problem that has remained unsolved in the strong-coupling limit is 
that of the so-called mode conversion: the semiclassical description breaks down in those points (or 
subspaces) of the classical phase space where the spin-orbit interaction locally becomes zero. 

In the present paper we shall apply the above methods to various two- and three-dimensional 
model systems, test their ability to reproduce the coarse-grained quantum- mechanical level densities 
of these systems, and explore their limitations. We shall also discuss the mode conversion problem 
that arises in the strong-coupling limit. Some preliminary results of our investigations have been 



presented in |21, 22 



Our paper is organized as follows. After a short reminder of semiclassical trace formulae for 
coarse-grained quantum systems in sec ^, we review in sec ^ the approaches of 17, [l^, 20 which 
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we then apply in the following. In sec 4.1 we investigate the two-dimensional electron gas (2DEG) 
with a spin-orbit interaction of Rashba type in an external magnetic field. In sec [4.2| we add 
to this system a laterally confining anisotropic harmonic-oscillator potential as a model for an 
anisotropic semiconductor quantum dot. Section |5| is devoted to a three-dimensional harmonic 
oscillator potential with standard spin-orbit interaction of Thomas type, as in the shell model for 
light atomic nuclei (see, e.g., ||23||). In sec [6|, we finally discuss the mode conversion problem and 
present some supporting evidence for the diabatic spin-flip hypothesis proposed in fl9|| . 

2 Trace formulae with coarse-graining 

The primary object of the semiclassical trace formulae is the level density (or density of states) 

g(E)=J2S(E-E k ) (1) 

k 

of a system described quantum-mechanically by the stationary Schrodinger equation 

Hcf> k = E k ^ k . (2) 

We stay throughout in the non-relativistic limit and assume the energy spectrum {E k } to be discrete. 
Classically, the system is described by a Hamilton function H (r, p) = E and the equations of motion 
derived from it. The function H(r, p) may be understood as the phase-space symbol of the operator 
H in the limit H — > 0. 

The level density (|l|) can be written as a sum of a smooth part and an oscillating part: 

g(E) = g(E) + 6g(E) . (3) 



The smooth part g(E) is in the POT given by the contribution of all orbits with zero length |25|| . It 
is often easily evaluated by the (extended) Thomas- Fermi theory or by a numerical Strutinsky aver- 
aging [26 1 of the quantum spectrum (see Jl5| for the relation of all three methods). The oscillating 



part Sg(E) is semiclassically approximated by trace formulae of the form 

5g sc {E) = J2-A P o(E) cos (\s po (E) - |a po J . (4) 

The sum is here over all periodic orbits (po) of the classical system, including all repetitions of each 
primitive periodic orbit (ppo). S po (E) is the action integral and a po the so-called Maslov index of a 
periodic orbit. The amplitude A po (E) depends on the integrability and the continuous symmetries 
of the system. When all periodic orbits are isolated in phase space, the amplitude is given by [jl| 

Ao(E) = \ . Sp (5) 
7r/2 ^|det(M po -l)| 

where T ppo is the period of the primitive orbit and M po the stability matrix of the periodic orbit. 
Examples of amplitude factors for systems with continuous symmetries or for integrable systems 
may be found in the literature quoted in the introduction. 

The po sum in (^|) does not converge in most cases; it must in general be understood 
asymptotic series that is only semiconvergent. However, much practical use can be made of trace 
formulae if one does not attempt to obtain an exact energy spectrum (given by the poles of the 
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level density), but if one is interested only in the coarse-grained level density. For this, we define a 
smoothed quantum-mechanical level density by a convolution of (|l]) with a normalized Gaussian: 

^^E'-^'w. (•) 

Here 7 is a measure of the desired energy resolution. Applying the convolution to the right-hand side 
of (|3|) will, for small enough 7, not affect the smooth part g(E). The convolution of the oscillating 
part, applied to the semiclassical trace formula (Q) and evaluating the integration as usual in the 
stationary-phase approximation, leads to the the coarse-grained trace formula |7], |27j 

5g sc (E) = £ e-W*°WApoW cos Q-M^) " ■ (7) 

The only difference to (Bj) is the additional exponential factor which suppresses the contributions 
from orbits with longer periods. Due to this factor, the periodic orbit sum now converges for not 
too small values of 7. Our choice of the Gaussian function in @ is rather arbitrary; cf [^TJ] where 
the regularization of the trace formula is discussed in terms of a general smooth test function. In 
many physical systems, experimentally observable quantum oscillations could be well approximated 
through such coarse-grained trace formulae in terms of only a few short periodic orbits (see fi"5| , 



16 [ for examples). Balian and Bloch |q] have used a small imaginary part of the energy, which 



corresponds to using a Lorentzian smoothing function. 



3 Periodic orbit theory with spin degrees of freedom 

In our present study, we want to apply the POT to systems of fermions with spin s = 1/2, in which 
the spin degrees of freedom are involved through an explicit spin dependence of the Hamiltonian. 
We write it as 

£ = £ 1 + #1 , H = ?- + V(r) , (8) 

2m 

and assume the spin-dependent part to have the following general form of a spin-orbit interaction: 

ffi = kC(r,p)-(7. (9) 

Here a = (a x , o~ y , a z ) is the vector built of the three Pauli matrices and 1 is the unit 2x2 matrix 
acting in the spin space spanned by the Pauli spinors. The Planck constant h in @ comes from the 
spin operator s = ^ha. The constant k is such that the spin-orbit term has the correct dimension 
of an energy; it is composed of natural (or material) constants but does not contain h. The vector 
C(r,p), which is the phase-space symbol of the operator C(r, p), may be interpreted as an internal 
magnetic field with arbitrary dependence on the classical phase-space variables r, p. In the standard 
non-relativistic reduction of the Dirac equation, C(r, p) becomes 

C(r,p) = [W(r)xp], (10) 

with k = l/4m 2 c 2 ; V(r) is an external electrostatic potential. For the spherical Coulomb potential, 
(||) and (|l0|) yield the Thomas term (corrected by a factor 2) well-known in atomic physics. 

It is by no means trivial, now, how to define a classical Hamiltonian corresponding to (||), since 
there is no direct classical analogue of the spin. Whereas numerous attempts have been made over 
the last seven decades or so to describe the spin classically or semiclassically, only two approaches 
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have lent themselves to an inclusion of spin degrees of freedom in semiclassical trace formulae. These 
are the approaches developed by Littlejohn and Flynn |17L [T^| , with extensions by Frisk and Guhr 
[19 1, and of Bolte and Keppeler [ 20 1 . We refer to these original papers for all details, as well as for 
exhaustive references to the earlier literature. Here we shall briefly present the resulting formulae 
which will be applied and tested for various model systems in the following sections. 



Bolte and Keppeler [ 20 1 started from the Dirac Hamiltonian to derive a relativistic trace formula 
which, to our knowledge, has not yet been applied to physical systems. They also started from the 
non-relativistic Pauli equation for a charged particle with spin 1/2 in an external magnetic field 
B(r), for which we have 



1 

2m 



P 



5a W 



+ V(r) , Hi 



eh 
2mc 



B(r) 



(11) 



Using the same techniques as in the derivation of their relativistic trace formula, they discuss 
two limits for introducing the semiclassical approximation. The Zeeman term Hi in ( |TT| ) is not a 
spin-orbit interaction, but Bolte and Keppeler [20| argue that the extension of their methods to 
the more general form (^) is straight forward. We therefore present their approach below for the 
general spin-orbit Hamiltonian (0). 



3.1 Weak coupling limit 

In the "weak coupling" limit (WCL), the semiclassical approximation is systematically performed 
by the usual expansion in powers of h. Because of the explicit appearance of h in Hi, the limit h — > 
leads to the classical Hamiltonian H c [(r, p) = Hq(t, p) whose periodic orbits enter the trace formula. 
The spin degrees of freedom here are not coupled to the classical motion. Their contribution to the 
trace formula enters through the trace of a 2 x 2 matrix d(i) obeying the "spin transport equation" 

Ad(t) = -iK[C(r,p).<x]d(i), d(0) = l, (12) 
at 

to be evaluated along each periodic orbit r po (t), p po (t) found from Hq(t, p). This equation describes 
the spin precession about the instantaneous internal magnetic field C(r, p) along the periodic orbit. 
Using the solution of QI2D for each orbit, the trace formula is given, to leading order in h, by p0| ] 

5g sc (E) = Ao(E) tr d{T po ) cos (\s po {E) - |a po J , (13) 

where T po = dS po (E)/dE is the period of each (repeated) orbit. Since the periodic orbits are not 
affected by the spin motion, the only difference to the standard trace formula (Q) is the appearance 
of the spin modulation factor trd(T po ); all other ingredients are evaluated in the usual manner for 
the unperturbed periodic orbits of Hq. One may therefore consider this treatment as an adiabatic 
limit of fast spin motion and slow spatial motion r(i), but Bolte and Keppeler argue that this 
adiabatic assumption is not needed for the formula (|l3|) to be true to leading order in h. 

Through the factor trd(T po ) the contribution of a given periodic orbit depends on the overlap 
of the spin directions at the beginning and the end of its period. An orbit for which these two 
directions are identical has simply trd(T po ) = 2, whereas an orbit for which these directions are 
antiparallel does not contribute at all to the trace formula (|l3|). 

For self-retracing orbits, i.e., librations between two turning points in coordinate space, the 



spin precession described by (12) is reversed at each turning point, and hence the spin direction 
is brought back to its initial value after a full period. Such orbits therefore only acquire a trivial 
factor tr d(T po ) = 2 compared to the trace formula (Q). For systems which possess only self-retracing 
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periodic orbits (see, e.g., the examples in sees [| and ||), the formula (^) thus reduces to the trivial 
recipe of incorporating the spin by a simple degeneracy factor 2 in the level density, which cannot 
account for the spin-orbit interaction effects. 



3.2 Strong coupling limit 



To obtain the "strong coupling" limit (SCL), Bolte and Keppeler [2C] follow the philosophy of [18, IS] 
by absorbing the Planck constant h in (|ll|) into the Bohr magneton \x = eh/2mc, thus considering 
[i as a constant in the semiclassical limit h — > 0. Similarly, for the spin-orbit Hamiltonian @ one 
absorbs h into the constant R = fin. The fact that this corresponds to a double limit h — > and 
k — > oo, with R = Hk kept constant, justifies the name "strong coupling" limit. 

The symbol of the full Hamiltonian in phase space now remains a 2 x 2 matrix which after 
diagonalization leads to the two classical Hamiltonians 

ff±(r J p) = J ff (r,p)±«|C(r,p)|, (14) 

which can be considered as two adiabatic Hamiltonians with opposite spin polarizations. They 
create two classes of dynamics, whose periodic orbits must be superposed in the final trace formula. 



Such a trace formula has, however, not been derived explicitly so far. Littlejohn and Flynn [17] 



argued that a non-canonical transformation of the phase-space variables r, p would be necessary to 



calculate the amplitudes. Frisk and Guhr [19] surmised, based upon Fourier transforms of quantum 
spectra, that this is not necessary, provided that the actions S^ a of the periodic orbits generated by 
the Hamiltonians H± be corrected by some phases accumulated along the periodic orbits: 



\ S% - ~ S± + A* ± , A* ± = + A™) tit . (15) 



The phase velocities A±, A^ N , which have been called the "Berry" and the "no-name" terms [17, 15], 
arise as first-order h corrections in the semiclassical expansion of the symbol of the Hamiltonian 
matrix. Bolte and Keppeler |2^] have used their techniques to give this prescription a rigorous 



justification. For the Hamiltonians (|i~4] ) with (|i~0|), the above phase velocities can be calculated 
most easily in terms of the polar angles 9, cj) defining the unit vector of the instantaneous direction 
of C(r,p), i.e., ec = (cos (p sin 9, sin cf> sin 9, cos 9), and of the Hesse matrix of the potential, V^- = 
d 2 V(r)/dridrj = x,y,z), evaluated along the periodic orbits, and are given by |l9| 



= T \ (1 - cos 9) <j> , \ N ± N = -^e T c V"e c . (16) 

Clearly, in the SCL the spin affects the classical dynamics, albeit only in an adiabatic, polarized 
fashion. Moreover, there is a serious limitation to the procedure outlined above. Whenever C = 
at a given point in (or in a subspace of) phase space, the two Hamiltonians H± become degenerate 
and singularities arise, both in the classical equations of motion and in the calculation of the phase 
corrections ( |i~6P and the stabilities of the periodic orbits. Such points are called "mode conversion" 
(MC) points. A similar situation occurs in the chemistry of molecular reactions when two or more 
adiabatic surfaces intersect. The MC poses a difficult problem in semiclassical physics and chemistry, 
that has not been satisfactorily solved so far for systems with more than one spatial dimension (see 



[28 1 for a discussion of the MC in one dimension). 

For self-retracing periodic orbits, all components of the momentum p are zero at the turning 
points. Hence, a spin-orbit interaction of the standard Thomas type (|i~0| ) will in the SCL lead 
to MC at the turning points. In this case, both the WCL and the SCL break down, and an 
improved treatment becomes necessary to include such orbits in a semiclassical trace formula. A 
new approach that is free of the MC problem has just been proposed [^] ; its results will be presented 
in a forthcoming paper. 
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4 Two-dimensional electron systems with Rashba term 



In this section we shall investigate a two-dimensional electron gas (2DEG) with a spin-orbit inter- 



action of the Rashba type [ 30 1 . We will include also an external magnetic field and an external 
potential V(r) = V(r x ,r y ) which causes a lateral confinement of the 2DEG (see, e.g., |H]]). The 
Rashba term can be written in the form 

/ -«)Pv \ 

H 1 = HkC(t,p)'&, C(r,p)= (v' z )p x . (17) 

\ p y dV(r)/dr x - p x dV(r)/dr y ) 

Here (v' z ) is the mean gradient of the electrostatic potential in the z direction that confines the 
electron gas to the (x, y) plane, so that z = p z = 0, and the constant k depends on the band 



structure of the semiconductor in which the 2DEG is confined [31]. Note that the Rashba term (|l 
is of the standard form (|l0|). 



We shall in sec 4.1 study the case of the free 2DEG for V(r) = in an external perpendicular 
homogeneous magnetic field, for which a quantum-mechanically exact trace formula is easily derived. 
We shall see that the WCL approach yields an analytical semiclassical trace formula which is exact 
only to leading order in fin and in the limit k — > 0. The SCL approach, however, for which we also 
obtain an analytical result, is demonstrated to fail for k — > 0, but to include correctly the higher- 



order terms in Hk, and to become exact in the strong-coupling limit. In sec 4.2, we add an anisotropic 
harmonic confinement potential V(r) to this system, where the WCL can be applied successfully 
in numerical calculations. When the external magnetic field is switched off, the remaining system 
has only self-retracing orbits for which both the WCL and the SCL fail. The MC problem arising 
in the SCL for this system will be discussed in section ^. 

4.1 Free 2DEG with Rashba term in an external magnetic field 



We first discuss the free 2DEG with the Rashba spin-orbit interaction ( |17| ) in a homogeneous 
magnetic field B = Bqg z . The corresponding vector potential A is included by replacing the 
momentum operator in the usual way: p — > 7r = p — eA/c. The total Hamiltonian then reads 

7T 2 ^ /-«>%\ 

# = — l + ^C(r,7r)-<7, C(r,iv)= I (v' z )n x J , (18) 

where, using the symmetric gauge for A, the mechanical momentum is given by 

*=(l)- e -^(~rl). (19) 



Here m* is the effective mass of the electron and e its charge. We have omitted the spin contribution 
to the Zeeman term which could be trivially included by adding the magnetic field to C. 

The quantum- mechanical eigenvalues of the Hamiltonian (111) are known analytically p0|: 



Px \ 




(~ r v 


Py 


l-f 1 


r x 


J 


\ o 



E = h J j c /2, E± = Huj c (n± ^1/4 + H2nK 2 ) , n = 1,2,3,... (20) 

Here uj c = eBo /m*c is the cyclotron frequency, and k = k (v ' z ) <J m*/ uj c is a renormalized coupling 
constant which we have defined in such a way that the Planck constant h appears explicitly in all 
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our formulae Q The exact level density is then given by 



g(E) = 5{E -E ) + J2 l S ( E ~ E n) + S ( E -En)}- 



(21) 



n=l 



Using Poisson summation (see, e.g., [15], sec 3.2.2), this result can be identically transformed to 
an exact quantum-mechanical trace formula. The smooth part of (^Tj) is g(E) = 2/huj c , and the 
oscillating part becomes 



Sg(E) 



oo 

£ 



oo 

< \ COS 

k=l 



y/l/A + 2EK 2 /uJ c + h 2 K i 

E 



k2ir{ — + Hk 2 ± a/1/4 + 2Ek 2 Iuj c + 
nuj r 



(22) 



We will now analyze the system semiclassically, using both the WCL and the SCL approach. In 
the weak-coupling limit, we first need the trace formula for the unperturbed system without spin- 
orbit coupling, corresponding to Hq = -k 2 /2m*. This is the quantized Landau level system, whose 
exact trace formula is that of a one-dimensional harmonical oscillator with the cyclotron frequency 
uj c and reads |l5| (without spin degeneracy factor) 



Sg 



(*=o) 



c fe=i v 



2ttE 

hjjj c 



(23) 



For the Rashba term in (|l8|), the spin transport equation (^) can be solved analytically [21|, and 
the spin modulation factor becomes 



trd(kT po ) = (-l) fc 2cos k 2ir v / l/4 + 2Ek 2 / 



lUJr 



With (|13|), the complete semiclassical trace formula in the WCL can therefore be written as 



Sg 



WCL 



< £ )^EE 



COS 



± k=l 



(El 
— ± J1/A + 2Ek 2 / 



LO, 



(24) 



(25) 



This result is not the same as the exact quantum- mechanical one (22), but it contains the correct 
terms of leading order in h, in accordance with the derivation of Bolte and Keppeler [20|, and the 
correct leading-order term in k 2 . The missing terms would come about by going to higher orders 
in the semiclassical h expansion. Note also that (25) becomes exact in the limit k — » 0. 

Recently, Keppeler and Winkler [32] have analyzed the anomalous magnetoresistance oscillations 
of a quasi-2DEG in GaAs semiconductors, employing two kinds of spin-orbit interactions one of 
which was of the Rashba type (^). They applied the WCL trace formula (|l3j) and obtained good 
agreement with quantum-mechanical results. As the spin-orbit interaction in GaAs is rather weak, 
we assume that their results were not sensitive to the missing higher-order semiclassical terms, 
which explains their good agreement. 

It is very instructive now to compare the above result with that of an analysis using the strong- 
coupling limit. The SCL Hamiltonians (|l4|) become H±(r,p) = Ho(r,p) ± K^2uj c Ho(r, p) with 
k = Kk. It is easy to see that Hq(t,p) = Eq is a constant of motion. The equations of motion 
derived from H± therefore become linear, representing one-dimensional harmonic oscillators like for 



1 In the literature on the Rashba term, the constant a = h 2 n(v' z ) is frequently used, see |30], 
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the Landau orbits of the unperturbed system H . Instead of the cyclotron frequencies tu c they have, 
however, the modified eigenfrequencies 



u± = ll> c ( 1 ± Ky/uc/ZEo) ■ (26) 



Since C(r, p) does not change its sign along the modified Landau orbits, the system does not suffer 
from the mode conversion problem and the SCL can be safely applied. The action integrals of 
the primitive orbits are simply found to be S± = 2ttEq/uj c , like for the unperturbed system, but 
expressing them in terms of the conserved total energy E = Eq ± H\/2uj c Eq one finds 

S±(E) =2n(J^ + K 2 T\j2E% 2 /uj c + K^ . (27) 

The phase velocities ( |T6| ) are easily found to be A± = =F0/2 = ±uj±/2 and \± N = 0, so that 
the overall phase correction (^) becomes AQ± = ±7r for each repetition of the primitive orbits. 
Inserting these results into the trace formula of the one-dimensional harmonic oscillator, i.e., eq 



(23) with ijj c replaced by w±, using T± = 2tt/uj± = dS±/dE, and summing over both orbit types, 



we obtain the semiclassical trace formula in the SCL (exhibiting again the h contained in k) 
sc V ' hwc^y J2EK 2 /u c + h 2 K A ) 



± 

<x> 



\ ( E X 
x^cos k2A — + Hk 2 ± v / 2Ek 2 /lu c + h 2 K 4 J 
k=i L V c /. 



(28) 



It is interesting to note that hereby the Berry term, yielding the phase correction kA<&± = ± kir, 
cancels the alternating sign (— l) k in (^). We note that the result (p8|) would be exactly identical 
to the quantum-mechanical result (|22|), were it not for the missing term 1/4 under the roots. We 
see, therefore, that the SCL result will fail in the limit k — ► 0, since the alternating sign (— l) k arises 
precisely from that missing term 1/4 in the actions. On the other hand, the SCL trace formula 
(p8|) does correctly include the higher-order terms in Hk 2 , both in actions and amplitudes, becoming 
exact in the limit of a large spin-orbit coupling parameter k, as could be hoped. Note that (|28|) 
becomes exact also in the limit of large energy E. That the SCL trace formula fails in the limit 
K — > is not surprising because of the non-analytic behaviour of the Hamiltonians (0), as already 



pointed out in [17]. 



4.2 A quantum dot with external magnetic field 

We now add a lateral confining potential V{r x ,r y ) to the previous system. This is a simple model 
for a two-dimensional quantum dot which nowadays can easily be manufactured in semiconduc- 
tor heterostructures. We choose the confining potential to be harmonic, so that the Hamiltonian 
becomes 

H = ^1 + ^ {u>l rl + u? y rl) 1 + Hk C(r, tt) • &. (29) 

The Rashba term C now acquires also a z component like in ([H]) and reads: 

-«>% 

C(r.7rj-[ _ (v' z )n x _ |. (30) 

m*ulr x 7t y - m*uj 2 r y TT x 



In the case where the two oscillator frequencies tu x and ui y are identical, the total system has axial 
symmetry and is integrable even including the spin-orbit term, with the eigenvalues of the total 
angular momentum in z direction, J z = L z + s z , being constants of the motion. An exact trace 
formula can then be found by the EBK quantization following the methods of |l^]. We will not 
discuss the integrable system here and refer the interested reader to |21|| . A less trivial situation 
arises when the frequencies uj x and lo v are different and the system with spin-orbit coupling is no 
longer integrable. 

Since a realistic spin-orbit coupling in most semiconductors is weak, wee shall here only use 
the WCL to derive a semiclassical trace formula for the Hamiltonian (^). The system without 
spin-orbit coupling is biquadratic in the space and momentum variables and can be transformed to 
become separable in its normal modes. The normal-mode frequencies are 



2 | 2 
J x + u y 



^/(w2+ W 2 +w 2)2_ 4aJ 2 w 2 



1/2 



(31) 



The exact eigenenergies are thus given in terms of two oscillator quantum numbers n + and n_ 



hw + (n + + 1/2) + huj-(n- + 1/2) n± = 0, 1, 2, 



(32) 



The semiclassical trace formula for such a system is known [ 33 1 and quantum-mechanically exact 



^ 1 " 0, < E > = £+>£(-* 



± =t fc=l 



sin [kn(uJzf/uj±)] 



sin ( k 



2irE\ 

hu>± J 



(33) 



Note that this formula is only useful when w + /u;_ is irrational. For rational frequency ratios, careful 



limits must be taken to cancel all singularities, see [33 1 for details. However, for a finite external 
field i?o 7^ 0, the ratio u; + /u;_ can always be made irrational by an infinitesimal change of the field 
strength, so that equation (|33|) is adequate for all practical purposes. The semiclassical origin of 
this trace formula is given by the existence of only two isolated rotating orbits with frequencies u;+ 
and uj-, whose shapes in coordinate space are ellipses. Each orbit contributes one of the above two 
sums; k is the repetition number of the primitive orbits (which have k = 1). 

We next have to calculate the spin modulation factors by solving the equation (12) along the 
unperturbed elliptic orbits. This could only be done numerically. It is, however, sufficient to 
calculate the modulation factors for the primitive orbits only. Using the property tr d(/cT po ) = 
trd fc (T po ), the final trace formula in the WCL is then given by: 



£ 



1 



oo 

y 



trd fc (Ti 



sin [kir(u> T /uJ±)] 



sin [k 



2itE 



(34) 



where T± = 2tt/u± are the periods of the unperturbed primitive orbits. 

In fig H| we compare results for the oscillating parts 5g{E) of the coarse-grained level density 
for 7 = 0.3 Hujq; all energies are in units of Hujq. The deformation of the confinement potential was 
fixed by uj x = ujq and u y = 1.23 wo, and the cyclotron frequency was chosen to be to c = 0.2 ojq. 
In our numerical calculations we have set H = uq = m* = e = c = {v' z } = 1. In these units, the 
spin-orbit coupling parameter was chosen to be k = 0.1. The heavy solid lines in the upper three 
panels represent the full quantum-mechanical result obtained from an exact diagonalization of the 
Hamiltonian (B^) in the basis of Ho. In the top panel, the semiclassical trace formula (|33|) without 
spin-orbit interaction is shown (only k = 1 and k = 2 contribute visibly). It clearly demonstrates 
that the effect of the spin-orbit interaction on the level density, even at this resolution, is dramatic. 
In the next two panels, the spin-orbit interaction has been included into the semiclassical WCL trace 
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formula (34), using the numerical spin modulation factors. We see that the agreement is improved 
radically, especially if the second repetitions (k = 2) are added. The difference between quantum 
mechanics (QM) and semiclassics (SC) can clearly be seen only in the close-up (second panel), which 
selects the energy region 11 < E/(?UjJq) % 21, where the disagreement actually is worst. The bottom 
panel shows the energy dependence of the two spin modulation factors tr d(T+) and tr d(T_) of the 
primitive orbits. Clearly, the strong long-range modulation in the amplitude of 5g(E) is the result 
of the spin-orbit interaction; it is correctly reproduced in the WCL approach through the inclusion 
of the spin modulation factors. 




^3 



SC (k=l) 
SC (k=l,2) 
QM 





Figure 1: Upper three panels: Coarse-grained (with 7 = 0.3 Hujo) oscillating part 6g(E) 
of level density of the 2-dimensional quantum dot with Rashba term (energy units: Hujo). 
Heavy solid lines: full quantum results (QM) including Rashba term. Dotted line in the 
top panel: trace formula ( |33|) for k = 0. Solid and dashed lines in 2nd and 3rd panel: 
semiclassical trace formula (|34|) (SC) with only first (k = 1) and up to second harmonics 
(k = 2) included. Lowest panel: spin modulation factors trd(T+) and trd(T_). 



This concludes the discussion of the system with magnetic field Bq ^ 0. We notice at this point 
that the case Bq = with irrational frequency ratio uj x /uj y is not accessible in the WCL. This 
follows from the fact that in the system without spin-orbit coupling, the only periodic orbits are 



the self-retracing librations along the principal axes. As already discussed at the end of sec 3.1, the 



WCL fails here in that it gives only the trivial modulation factor trd(T po ) = 2 for both orbits. On 
the other hand, the SCL suffers from the MC problem. We will discuss this problem explicitly in 
sec H|, where we return to the above system with Bq = 0. 
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5 Three-dimensional harmonic oscillator with standard 
spin-orbit interaction 

We now discuss a three-dimensional system with a spin-orbit interaction of the Thomas type, as it 
is well-known from non-relativistic atomic and nuclear physics. In order to be able to perform as 
many calculations as possible analytically, we choose again a harmonic-oscillator potential for V(r). 
This potential is not only the prototype for any system oscillating harmonically around its ground 
state, but it has actually been used in nuclear physics as a realistic shell mo delQfor light nuclei || ; 
provided that the spin-orbit interaction was included with the correct sign 
We thus start from the following Hamiltonian 



H = H 1 + HkC(t,-p)-&, C(r,p) = Vy(r)xp, (35) 



with 



H = ^ + V(r), V(r)= £ ^ 2 rf. (36) 

i=x,y,z 

Here r = (r x ,r y ,r z ) and p = (p x ,Py,Pz) are three-dimensional vectors. We express the three 
oscillator frequencies in terms of two deformation parameters a, /3: 

lv x =uj , Lo y = (1 + a) uj , Lo z = (1 + a) p u) Q . (37) 

and use Hioq as energy units. For a = 0, j3 = 1 the system has spherical symmetry, for (5 = 1 and 
a 7^ it has only axial symmetry, k will be measured in units of (Huq) . 



We shall first (sec 5.1) briefly discuss the quantum-mechanical spectrum of the system, and 
then (sec |5.2| ) investigate it in more detail by the semiclassical methods. The most interesting case 
is that where the oscillator frequencies are mutually irrational, so that the unperturbed classical 
Hamiltonian Hq has only self-retracing periodic orbits. In this case the WCL cannot handle the 
spin-orbit coupling, and we must resort to the SCL. As we will show, the leading periodic orbits 
with shortest periods in this case do not undergo mode conversion. We therefore use this system 
for a representative case study, for which a trace formula can be successfully derived (cf |22|, |2l[ |) 
within the SCL. 



5.1 Quantum- mechanical spectrum 

In general, the system ( p5| ) is not integrable. There are, however, two well-known integrable cases 
for which the quantum spectrum is analytically known: the separable system (^) without coupling 
(k = 0), and the spherical system (a = 0) including coupling. The unperturbed harmonic oscillator 
has the spectrum 

En x n y n z = ^ ^ ^ + l l 2 )' ^ = 0, 1, 2, . . . (38) 

i=x,y, z 

The spherical system (a = 0) with coupling for which C(r, p) = ojq L, has the spectrum 

E nlj =hnj {2n + l + Z/2)+K{hw ) 2 x{_^ +i) ^ 3 . = { + ' (39) 

2 The angular-momentum dependent L 2 term included in the Nilsson model ]2g| ] is of minor importance 
in light nuclei; it is left out here to simplify our investigation 
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where n,l = 0, 1, 2, . . . and j = Z± 1/2 is the total angular momentum. The spin-orbit term is only 
to be included for / > 0. Each level E n ij has the usual angular momentum degeneracy (2j + 1) 
which equals 2 for I = 0. Note that in nuclear physics, k is negative [p3, 24]. 



The non-integrable cases require numerical methods for determining the energy spectrum. Here 
we used the diagonalization in the basis i Uy \n x , n y , n z , s z ) of the unperturbed Hamiltonian Hq ( p6| ) 
with the eigenenergies (|3~8|), where \s z ) with s z = ±1 are the spin eigenstates of a z . The inclusion 
of the phase i ny leads to real matrix elements; furthermore, the conservation of the signatures 
(— X s j n *+ n y+ n ' a nd (—l) nx+ny s z allows one to separate the Hamiltonian matrix into smaller uncoup- 
led blocks (see Q for details). 



5.2 Semiclassical analysis 
5.2.1 Smooth level density 

When one wants to compare results of semiclassical trace formulae with quantum-mechanical level 
densities, one has to subtract from the latter the smooth part g(E) (see sec S). For the Hamiltonian 



(35), g(E) can be calculated analytically within the extended Thomas-Fermi (TF) method, which 



has been done already long ago [35]. The result, as an expansion both in ft and powers of k, reads 
E 2 



9(E) = ^ {l + t?« l { l 4+J} + u 2 z ) + 0(h* K *)\ 

ft? uj T uj„ujy I " J 



' U X UyUJ Z 

2E 



+ 77T2 i"Wy + "Wz + + 0(K 5 K*)\ 

6 ft z u t w„uj z I J 



' U X UJyLV Z 

(uj 2 x + + u? z ) f | ^ (aj + oj + u; 2 z ) 2 + 2 frjoj + Jfal + ujoj) | ^,,4)] _ (4Q) 
12hu x u y u z \ (to 2 + uj 2 + uj 2 ) )' 

In the literature, the smooth part is often assumed to be given by the TF model. This leads, 
however, only to the leading term proportional to E 2 . For an accurate determination of g(E), the 
leading ft and ft 2 corrections relative to the TF term may not be neglected. 

5.2.2 Trace formulae for the integrable cases 

The exact spectra of the integrable cases offer again the possibility to derive trace formulae that are 
exact in all orders of ft. For the unperturbed harmonic oscillators, these are known |l5|, |33|] and need 
not be repeated here. For the spherical case with spin-orbit interaction, the methods of [|l5|, 33, 35] 
lead to the following result: 

x (j* E 1 1 . f u ET± 7T L 

9[ ' ~ (ftw ) 2 ^^ i (l±Kftw ) 2 sm[2k7r/(l±Kftio )} Sm \ * * 

1 ^y, (±1 + 2^q) 1 zh J h ET ± l ka+ 

ftw 2(1 ±Kftuj ) 2 sin[2A;7r/(l ± kHojq)] \ ft 2 

1 v-^ v~~^ 1 cos [2/c7r/ (1 i KftidQ)\ ( ET± ir 

| X x COS I K K(T l 

ftuj z ^f^ i (l± k?uv ) 2 sin 2 [2kn/(l ± kHu )] \ ft 2 



+ 



ftw ^ 2sin 2 (A:7r^o) V & 2 / 
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where the three periods T± and Tq are given by 



ujq(1 ± kHluq) ' UJq 

and the phases 



To = * (42) 



±2 

cr± = — — - — , cr = -Anfrwo (43) 
1 ± kHluq 



play the role of non-integer Maslov integers. When added to the smooth part (p0|), eq (^l|) repro- 



duces the exact quantum spectrum (3£). This trace formula thus serves us as a test limit of the 
semiclassical results derived below in the non-integrable deformed cases. In To we recognize the 
period of the classical orbits of the unperturbed Hamiltonian; the shifted periods T± have to be 
explained by the periodic orbits of the perturbed system. 

In the limit k — > 0, the sum of the smooth term ( (40| ) with uj x =u> y =u! z =u>o and the oscillating term 
( f4l"D yields the exact trace formula of the isotropic three-dimensional harmonic oscillator |I5| , |33| 
(which here includes the spin degeneracy factor 2) 



9(E) ' 



E 2 - \(f^o) 2 



jl + 2 ^(-l) fc cos (k 
^ k=i ^ 



2-kE 



(44) 



5.2.3 Fourier transforms 

Since we have no explicit quantum spectra of the perturbed system and therefore cannot derive an 



exact trace formula, we resort to the method of Fourier transforms of the quantum spectrum [36] 
in order to extract information on the periods of the system. The Hamiltonian (^) has the scaling 
property 

H (rfr, 7/p) = 7] 2 H{r, p) . (45) 

We see below that this scaling property holds also in the classical limit if the SCL is used. As 
a consequence, the energy dependence of the classical dynamics and thus of the periodic orbits is 
simply given by a scaling, and their actions go like S po (E) = T po E, whereby the periods T po = 2tt/uj po 
are energy independent (but depend on k). Therefore, the peaks in the Fourier transforms of 5g(E) 
with respect to the variable E will give us directly the periods T po in the time domain, whereby the 
peak heights are given by the semiclassical amplitudes A po and their signs give information on the 
relative Maslov indices. 

In fig ^ we present a series of Fourier transforms of 5g(E) obtained from the numerically diago- 
nalized quantum spectra with a coarse-graining parameter 7 = 0.5 Hujq. Shown are here the squares 
of the Fourier amplitudes in the time domain, plotted for different spin-orbit coupling strengths n. 
A slightly anisotropic ratio of frequencies u x = ujq, oj y = 1.1215wo, uj z = 1.2528 a>o was chosen. For 
k = 0, the system then has only the three isolated librating orbits along the principal axes. Indeed, 
wee see at n = the three dominant peaks with the corresponding primitive periods Tj = 27r/a;j 
(i = x, y, z). Their second harmonics (k = 2) are also resolved; however, due to their larger periods 
they are of smaller amplitude. For k > this simple peak structure is split and ends in a completely 
different spectrum at k = 0.2 (Hloq)' 1 . The dotted lines are the predictions from the semiclassical 
SCL analysis given in the next subsection; the plots on the right-hand side indicate the shapes of 
the periodic orbits for increasing k. 
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5.2.4 Semiclassical treatment in the SCL approach 



We now analyze the system semiclassically in the SCL. According to sec 3.2, the classical Hamilto- 
nians to be used are 



H ± (r,p) 



ip 2 + y(r)±R|C(r,p) 



(46) 



with C given in (p5[) and k = Hk as discussed in sec 3.2. The spin-orbit term destroys the integrability 



of the harmonic oscillator, but the scaling property (45) is still fulfilled. Therefore the above Fourier 
spectra should give us the correct periods of the periodic orbits defined by the Hamiltonians 



k=0.20 



k=0.16 




k=0.08 



k=0.04 



Figure 2: Fourier spectra of quantum-mechanical level density 5g(E) (coarse- 
grained with 7 = 0.5 Itwq) °f the 3-dimensional harmonic oscillator with defor- 
mation a = 0.1212, P = 2, and various spin-orbit strengths k in units (htoo)^ 1 . 
T is in units of uJq 1 . For the dotted lines and the inserts on the rhs, see fig ||. 



a) Periodic orbits 



The equations of motion for the Hamiltonians (46) become 

Pi ± e ijk R \C\~ x (Cjulr k - C k Ujrj) , 
-u^n ± e ijk R IC]" 1 (CjU^p k - C k u%pj) . 



Pi 



i,j,k G {x,y,z} 



(47) 
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This is a non-linear coupled system of six equations, and the search for periodic orbits is not easy. We 
have determined them numerically by a Newton-Raphson iteration employing the stability matrix 
[ 21 1 . Special care must be taken at the MC points where C(r,p) = and hence the equations (47) 
are ill defined. In general, this leads to discontinuities in the shapes of the periodic orbits, due to 
which their stabilities cannot be defined. We shall return to the MC problem in sec |(| It turns out 
that there exist periodic orbits which are free of MC, i.e., for which C(r, p) never becomes zero. The 
existence of some particularly simple orbits follows from the fact that the three planes r k = pj, = 
{k = x, y, z) in phase space are invariant under the Hamiltonian flow. The coupled equations of 
motion for the class of two-dimensional orbits lying in the planes are 



h = Pi =F eijkKsig^C^ujj rj , 

Pi = -ujfn =F e ijk R sign(C fc ) u>f pj , i, j, k G {x,y,z} 



(48) 



where i and j refer to the in-plane variables and k to the normal of each plane. Assuming that 
there exist solutions with Ck / 0, we can put sign(Cfc) = 1, since to each such orbit, there exists a 
time-reversed partner which belongs to the opposite of the two Hamiltonians H±. Hence we obtain 
the system of equations 



n 

Pi 



Pi 



-ufri =F eijkKujfpj 



i,3,k G {x,y,z} 



(49) 



which now are strictly linear and can be solved by finding the normal modes. The solutions are two 
periodic orbits of ellipse form in each invariant plane For the two Hamiltonians H± this gives 

altogether twelve planar periodic orbits which come in doubly degenerate pairs. We therefore find 
six different orbits (denoted by ^f- in fig || below) with the frequencies 



^'7 



UJ-; 



1/2 



(50) 



The rhs of fig |3| shows these frequencies versus the spin-orbit coupling parameter R. On the lhs 
we illustrate the periodic orbits for R / (ellipses) and the unperturbed libration orbits for R = 0. 
The periods = lixju)^ fit perfectly the positions of the most pronounced peaks in the Fourier 
spectra of fig |] (cf the dotted lines there), when the appropriate deformations uji, ujj are chosen. 
Some of the minor peaks may be attributed to non-planar orbits (see below). 



r k =0 



K > 



K = 




Figure 3: Right: the six frequencies uif- fl50| ) (units: ujq) for ij = 12, 23, 

and 31 versus R (units: Uq 1 )- Deformations as in fig Left: schematic 
plot of shapes of the elliptic orbits {R > 0) and the unperturbed librating 
orbits {R = 0) in the plane. 



In the spherical limit, the same procedure leads to a very simple analytical result. Due to the 
conserved angular momentum L = r x p, most periodic orbits are planar circles. In each plane, 
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the equations of motion are similar to (|49|), with the two eigenfrequencies uj^ and corresponding 
periods given by 

wo(l±«wo), T ± = _ . . (51) 



WO (1 ± K^o) ' 

The periods T ± are exactly equal to the two periods T± in ( ^) that appear in the exact trace formula 
( f4l|) of the spherical system.^] However, the unperturbed harmonic oscillator period To = 2it/ll!q 
that also appears in (|4l|) cannot be explained by the present solutions in the SCL. We surmise 
that it might be connected to the existence of straight-line librating orbits; these lie, however, on 
mode conversion surfaces and cannot be treated in the present approach. In our ongoing studies 
[p9f where the MC problem is avoided, we can, indeed, confirm the existence of periodic orbits with 
the period To. 

Besides the above harmonic planar solutions, the full non-linear system ( |47| ) leads also to non- 
planar three-dimensional periodic solutions with C(r, p) > (or < 0) for which mode conversion 
does not occur. Some of these numerically obtained orbits evaluated at R = 0.2 oJq 1 are shown in 
fig ||. We also give their periods T po and partial traces Aj which determine their stabilities (see 
below). The orbit 7^ 2 has the period T po = 10.08 Uq which seems to be supported numerically by 
a small peak seen in the uppermost Fourier spectrum of fig [2|. The periods of the other orbits from 
fig U could not clearly be identified in the Fourier spectra; some of these orbits are too unstable and 
some of the periods lie too close to those of the leading planar orbits. 
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7 243,B 



T po =15.35 oj 

Ai=-0.64 

A 2 =0.50 



-0.75 0.0 0.75-1.5 0.0 1.5 

T po =15.34^ 1 
A-i=-1.38 
A 2 =4.52 



Figure 4: Shapes of four 3-dimensional non-planar orbits found in the harmonic oscillator 
with deformations as in fig ^ and spin-orbit interaction R = 0.2 uJq 1 , projected onto the 
three spatial planes. (See text for the periods T po and stability traces Ai, Ag.) 



'Note that in the present SCL approach, the constant n includes a factor h, see sec 3.2 



16 



Stability amplitudes and trace formula for isolated orbits 

The amplitude A po of a periodic orbit in the trace formula is strongly dependent on its stability. 

For isolated orbits in a two-dimensional system, the factor | det(M po — 1)| in (|5|) equals |2 — trM po | = 

|2 — (Ai + A2) I = |2 — (Ai + 1/Ai ) I , where Aj {i = 1, 2) are the eigenvalues of the stability matrix M po , 

and thus the quantity trM po contains all information about the stability of an orbit. For a system 
in d > 2 dimensions, we can write 

d-l 



det(M po -l)| = JJ|Ai-2| 



(52) 



where the "partial traces" Aj are the sums of pairs Aj, 1/Aj of mutually inverse eigenvalues of the 
(2d— 2) dimensional stability matrix M po : Aj = Aj + l/Aj (i = 1, 2, . . . , d— 1). An orbit is stable when 
|Aj| < 2 for all i, unstable when |Aj| > 2 for all i, and mixed stable (or loxodromic) in all other cases. 
For the latter cases, the stability depends on the phase-space direction of a perturbation. Whenever 
Aj = +2 for any partial trace, a bifurcation occurs and the stability denominator (52) becomes zero. 
In such a situation one has to resort to uniform approximations [12, |ll|] in order to obtain finite 
semiclassical amplitudes. Non-isolated periodic orbits with Aj = 2 occur in degenerate families 
for systems with continuous symmetries and are characteristic of integrable systems; for these, the 
amplitudes must be obtained differently |], [?], ||, || . The symmetry breaking away from integrability 
can also be handled perturbatively JTQ] or with suitable uniform approximations [11, 14]. 

Most of the orbits that we have found, both planar and non-planar, undergo bifurcations when 
the spin-orbit parameter R or the deformation parameters a, f3 are varied. A typical scenario is 
illustrated in fig ||, where a is varied at fixed R = 0.1 uJq 1 and (3 = 2. Shown are the partial traces Aj 
of the involved orbits. One of the planar ellipse orbits (7^) lying in the (x, y) plane undergoes an 
isochronous pitchfork bifurcation at a = 0.3977. The new-born pair of orbits (7^) are degenerate 
with respect to a reflection at the (x, y) plane. They are non-planar warped ellipses which rotate 
out of the (x, y) plane when a is increased, and then approach the (x, z) plane. Through an inverse 
pitchfork bifurcation at a = 0.4450, they finally merge with another planar ellipse orbit (7^) lying 
in the (x, z) plane. Near a ~ 0.425, the orbit 7j~j f suffers from two more bifurcations (the other 
orbits involved thereby are not shown). 



X 




0.375 



0.4 
a=0.3977 



0.45 
a=0.4450 



Figure 5: Bifurcation of a planar ellipse orbit under variation of the deformation param- 
eter a. Upper panels: projection of the orbits onto the (y, z) plane. Lower panel: partial 
traces Aj of the involved orbits versus a (see text for details). 
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This example shows that the classical dynamics of the Hamiltonians H± is mixed and quite 
complicated due to the unavoidable bifurcations. In principle, isolated bifurcations can be handled 
using the by now well-known uniform approximations [12, 13]. These fail, however, when two bifur- 
cations lie so close that the difference between the corresponding actions S po becomes comparable 
to or less than h. It would lead outside the scope of the present study to attempt to regulate the 
Gutzwiller amplitudes by uniform approximations. In fig ^ we show by crosses the critical values 
of the frequencies uj v and uj z (in units of u x = ujq) versus R, at which bifurcations of the planar 
orbits - occur for fixed values of the deformation parameter (3 = 2 (left) and = 3 (right). 



The other deformation parameter a is given via (37). In the deformation regions below the dashed 
lines, where no bifurcations occur, the semiclassical amplitudes can be used without further uniform 
approximation . 
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Figure 6: Bifurcations of the planar periodic orbits 7^ in the 3-dimensional harmonic oscillator 
with spin-orbit interaction, given by the SCL Hamiltonians H± in (]16|). The crosses give the critical 
values of the frequencies u) y and lo z versus R for (3 = 2 (left) and (3 = 3 (right). No bifurcations 
were found in the regions below the dashed lines. 



For the contributions of all the isolated periodic orbits away from the bifurcations, we therefore 
use the trace formula 



po 



-( 7 T po /2ft) 2 



-ppo 



cos I \ S vo + A$p - Z cr 



po , 



(53) 



whereby the sum po explicitly includes all periodic orbits of both Hamiltonians H±. The Maslov 



indices a po were evaluated with the methods developed by Creagh et al [37], employing the recipes 
given in appendix D of fl5[| . The terms A$ po are the phase corrections (15). For the planar ellipse 



and \% N 



orbits lying in the planes, we find Aj 

In fig [7] we show the results obtained for the situation a 



po 



-Rul/2, so that A$ 
= 0.1212, /3 = 2,R = 0.1 uj q \ for which 
no close-lying bifurcations exist. The quantum-mechanical coarse-grained level density 5g(E) is 
shown by the solid lines (QM) and includes the spin-orbit interaction in both curves a) and b). The 
semiclassical results (SC) are shown by dashed lines; in a) without spin-orbit interaction, which 
again demonstrates that the latter dramatically changes the level density; in b) with spin-orbit 
interaction through the trace formula (|53"|). Only the six primitive planar orbits have been used. 
We see that this already leads to an excellent agreement with quantum mechanics, except at very 
low energies where semiclassics usually cannot be expected to work. The curve SC in the lowest 
panel c) shows the semiclassical result over a larger energy scale. 
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Similar results were obtained in the region of deformations and R values below the dashed lines 
in fig ^, where bifurcations do not occur. In all these cases, it turned out that the inclusion of the 
six primitive planar orbits was sufficient within the resolution given by the coarse-graining width 
7 = 0.5 frujQ. This result is in agreement with the Fourier analysis of sec 5.2.3 of the quantum 



spectra, where all dominant peaks correspond to the periods of these six orbits. 



a) 



b) 



c) 



""So 




5 10 15 20 25 30 35 40 45 




Figure 7: Coarse-grained level density Sg(E) of 3-dimensional harmonic oscillator with 
spin-orbit interaction R = 0.1 oJq 1 (other parameters as in fig |2|). Upper panels: solid 
lines (QM) give the quantum-mechanical results with spin-orbit interaction. Dashed lines 
(SC) give the semiclassical results according to (|53l), calculated a) without and b) with 
spin-orbit interaction. Bottom panel: c) same as SC in b), but over a larger energy region. 
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6 The problem of mode conversion 



In this section we want to discuss the mode conversion (MC) problem that arises in the strong- 



coupling limit (SCL) following [17, 18]. In particular, we will discuss an intuitive method, suggested 
by Frisk and Guhr to partially avoid the MC problem. This method can qualitatively explain 
some of the peaks observed in the Fourier spectra of the quantum-mechanical level densities 5g(E). 
It can, however, not be used to calculate the amplitudes A po required for the semiclassical trace 
formula. 

To this purpose, we return to the 2DEG with lateral harmonic confinement discussed in sec 
4.2, but without external magnetic field (Bq = 0). We shall again assume the oscillator frequencies 
oj x and uj y to be incommensurable (i.e., oJ x /^>y is irrational). Then, the system without spin-orbit 
interaction has only the isolated self-retracing librating orbits along the axes, and the weak-coupling 
limit (WCL) approach of j2(| cannot handle the spin-orbit interaction (except for a trivial spin factor 
2 in the trace formula). We therefore have to resort to the SCL approach. In order to simplify the 
discussion and focus on the important points, we set m* = 1 and ignore the diagonal elements of 
the spin-orbit interaction. We thus start from the Hamiltonian 

H = l(pl + f y + ulrl + u;y y )l + H K (^ , (54) 

which in the SCL leads to the classical Hamiltonians (R = Tin) 



1 



H± = - [p x +P y + u x r x + u y r y ) ±K^Jpl+p%. (55) 

Before taking the semiclassical limit, we first perform a Fourier analysis of the quantum spectrum 
of (|54j ) which is easily diagonalized in the unperturbed harmonic-oscillator basis. The Hamiltonian 
fl54l) does not possess the scaling property (|45|) of the three-dimensional system studied in sec [|. 



However, we can use the method of [36] by scaling the parameter k away. Dividing equation (|5 



by k 2 and introducing the scaled variables fj = rj/K, pi = pi/n, we obtain the scaled Hamiltonians 

H± = \ {fx +f y + w 2 x tI + wjfjj) ± yjp*+p* = E/R 2 = e , (56) 

so that the classical dynamics does not depend explicitly on R but is determined only by the scaled 
energy variables e. Therefore, a Fourier transform of the quantum spectra along the path in the 
(E, k) plane with constant E/k 2 leads in the time domain to the quasiperiods T^ Q = s^ a (e)/e of the 
periodic orbits of the Hamiltonians (|56|) , whereby s^ a are their scaled actions. 

In the lower part of fig || below we show the result of the Fourier transforms, taken at two different 
scaled energies. For e = 10 5 Hluq/R 2 (dashed line) only two peaks are seen. For e = 30frwo/R 2 (solid 
line), these are slightly shifted and remain the dominant peaks, whereas four additional small peaks 
appear (the second of these extra peak is almost absorbed in the left dominant peak). 



We now analyze the classical dynamics of the Hamiltonians (56). Like in the three-dimensional 



system analyzed in sec 5.2.4 , we can find orbits lying in the invariant subspaces of phase space with 
?i = Pi = for one of the degrees of freedom i (x or y). For the other degree, the one-dimensional 
equations of motion become 

fj = pi ± sign pi, 

Pi = -ufh- (57) 

On the axes pi = in phase space these equations are ill defined. It is still possible to solve the 
equations for pi / 0, which leads to portions of a circle in the (pi, fj) plane for each sign of pi. One 
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may then connect these partial trajectories to form periodic orbits whose shapes have, however, 
kinks. This is illustrated in the upper part a) of fig ||. Due to the kinks, these orbits are not 
differentiable and their stabilities cannot be defined. Their periods can, however, be calculated 
analytically and become 



T?f a (e) = — arccos 



1 - 2e 



lf^(e) = — -— arccos 



4-7T 



1 - 2e 



\ 1 + 2e / ' <jJi uJi \1 + 2e 

From the area enclosed by the orbits, we can also obtain the scaled actions: 



„adia 



(58) 



(59) 



We use the superscript "adia" because the Hamiltonians H± correspond to the adiabatic situation 
where the spin polarizations are fixed. 




Figure 8: Upper part: a) Periodic orbits found from the adiabatic 
Hamiltonians (56). Lower part: b) Diabatic periodic orbits found 
by enforcing spin flips H + *-* H_ at the mode conversion points. 



An alternative way to use the partial solutions found from ([57]) has been proposed by Frisk and 
Guhr [jXQ|| : instead of joining the two portions obtained for both signs of pi on one and the same 
of the Hamiltonians H±, one switches between the H±, enforcing a spin flip at the MC points: 
H + < — * H-. This corresponds to the transition from the adiabatic to a diabatic basis. The orbits 
thus obtained are continuous circles with continuous derivatives, as illustrated in the lower part b) 
of fig ^, and correspond to simple harmonic librations along the i axes. Their periods and actions 
are easily found to be 



j.dia _ 



S f a (e) 



1\ 2tt 
e + - — 
2 J ^ 



(60) 



The superscript "dia" indicates that we call these the diabatic orbits. Their periods are those 
of the Hamiltonian without spin-orbit coupling: Tf 1& = = 2ir/u)i. Although their shapes 
are continuous and differentiable, their stabilities still cannot be calculated, because the Hessian 
matrices of the H± in phase space are singular at the MC points. In the limit e — > oo, the pairs of 
adiabatic orbits merge into the diabatic orbits, and we get 



rriadia / 
1 i,± 



-idia 



<± a (e) 



ndia 



(61) 
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s(e)/e 




Figure 9: Lower part: Fourier spectra like in fig ^ of the quantum spectrum of the 
Hamiltonian ( |54|) with 7 = 0.2 Hojq, evaluated for two values of the scaled energy e 
(units: Hujo/k 2 ). Deformation: uj x = u>q, uj v = 1.22>ujq. Upper part: quasiperiods 



s(e)/e from ( |59| ) and ( pOj ) for adiabatic orbits (dashed lines) and diabatic orbits 
(solid lines), respectively, of Hamiltonians (Egf) versus scaled energy e. 



In the upper part of fig ||, we show the curves s a< ^ a (e)/e by the short and long dashed lines, 
and the curves sf ia (e)/e by the solid lines. We see that their values at e = 30hu>o/R 2 correspond 
exactly to the six peaks appearing in the corresponding Fourier spectrum. The two dominant peaks 
correspond to the diabatic orbits, and the four small peaks correspond to the adiabatic orbits. For 
e = 10 5 /iu;o/k 2 , the only two peaks correspond to the asymptotic values of sf ia (e)/e = Tj dia , in 
agreement with the limit (61). 

The evidence of diabatic orbits according to the above spin-flip hypothesis had already been 
observed by Frisk and Guhr Jl9[] . They did, however, not recognize any signatures in their Fourier 
spectra corresponding to periods of adiabatic orbits involved with MC points, such as we have found 
them in the four minor peaks of fig ||. Their conclusion was therefore that spin-flips always occur 
at the MC points. Our results seem to suggest that both kind of dynamics occur. The dominant 
Fourier peaks are, indeed, those corresponding to the diabatic orbits which undergo spin flips at the 
MC points. However, there must also exist a finite probability that the orbits stay on the adiabatic 
surfaces H± = E, leading to the smaller peaks positioned at the correct adiabatic quasiperiods 
s(e)/e. 

The semiclassical amplitudes required for the trace formula cannot be calculated for the present 
system, neither using diabatic nor adiabatic orbits. The fully polarized treatment of the spin 
variables used in the SCL approach is obviously not flexible enough to account for the full dynamics, 
although the Fourier analysis of the quantum spectra suggest that there is some partial truth to it. 
A more complete semiclassical description of the spin motion should allow for a balanced mixture 
of adiabatic and diabatic spin motion. 
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7 Summary and conclusions 



We have derived semiclassical trace formulae for several nonrelativistic two- and three-dimensional 
fermion systems with spin-orbit interactions of Rashba and Thomas type. We have thereby employed 
the weak-coupling limit (WCL) developed by Bolte and Keppeler [20 1 , and the strong-coupling limit 
(SCL) of Littlejohn and Flynn [17| with extensions and justifications of Jl^, |20| , 

In the WCL approach, the spatial motion of the particles is taken into account only using the 
periodic orbits of the system without spin-orbit interaction. The spin motion is included adiabat- 
ically via the trace of a spin transport matrix d(t) which describes the spin precession about the 
instantaneous magnetic field provided by the spin-orbit interaction. Nevertheless, we found that 
for a 2DEG in an external magnetic field with and without lateral anharmonic confinement, the 
WCL yields excellent results. In the free case, for which an exact trace formula can be derived, 
the semiclassical WCL reproduces the exact leading-order terms both in h and in the spin-orbit 
coupling constant k. In the laterally confined case, the gross-shell structure of the coarse-grained 
quantum level density was very accurately reproduced numerically. From our results it can be seen 
that in the limit of very large spin-orbit constants k, the missing higher-order terms may restrict 
the applicability of this method. A particular situation, where the WCL approach misses the effects 
of the spin-orbit interaction totally, is that where only self-retracing isolated periodic orbits exist, 
for which the trace of d(T po ) only yields a trivial spin degeneracy factor 2. 

We have studied the SCL approach for two systems possessing exclusively self-retracing isolated 
orbits for which the WCL approach fails, namely two- and three-dimensional harmonic oscillators 
with irrational frequency ratios. In the SCL approach the mode conversion (MC) problem, arising 
at points in phase space where the spin-orbit interaction locally is zero, imposes severe restrictions, 
since singularities in the equations of motion and the linear stability analysis of periodic orbits arise 
at the MC points. However, in the three-dimensional case, which provides a realistic shell model for 
light atomic nuclei, we found that the leading orbits with shortest periods are free of MC and lead 
to excellent results of the semiclassical trace formula for the coarse-grained level density, as long as 
bifurcations of these orbits are avoided. (The latter, when they are sufficiently separated in phase 
space, can be taken into account using well-developed uniform approximations and do, in principle, 
not affect the applicability of the SCL approach.) 

In the two-dimensional model of an anisotropic quantum dot, the MC problem could not be 
avoided. By a Fourier analysis of the quantum spectrum we have provided some support of the dia- 
batic spin-flip hypothesis put forward by Frisk and Guhr fll9|| . We have to extend their conclusions, 
though, in the sense that there is evidence for a mixture of both diabatic and adiabatic classical 
motion on the two spin-polarized energy surfaces H± = E, whereby the diabatic periods dominate 
the Fourier spectra. But even in the diabatic spin-flip limit, the semiclassical amplitudes of the 
trace formula cannot be calculated, and the MC problem therefore remains essentially unsolved. 
The connection between the existence of MC points and real spin-flip processes should therefore 
be taken with caution. In order to study the physical relevance of spin flips in the presence of a 
spin-orbit interaction, there is a definite need for a better analytical semiclassical treatment of the 
spin degrees of motion that is free of singularities and allows for a balanced mixture of adiabatic 
and diabatic spin motion. A new approach |^9[ in which the MC problem does not arise is presently 
being developed and will be presented in more detail in a forthcoming paper. 
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